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Abstract: We consider the problem of the Winterbottom's construction 
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microscopic validity within a 1 + 1-dimensional SOS type model. We then 
present the low temperature expansion of the wall tension leading to the 
Wcnzel's law for the wall tension and its corrections. Finally, for a fix rough- 
ness, we compare the influence of different geometries of the substrate on 
wetting properties. We show that there is an optimal geometry with a given 
roughness for a certain class of simple substrates. Our results are in agree- 
ment and explain recent numerical simulations. 
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1 Introduction 



Wetting phenomena have a long standing history starting with Young more 
than a century ago. His famous equation describes the behaviour of the 
contact angle 9 of a sessile liquid drop B in equilibrium with the vapor phase 
A on top of a substrate W: 

t ab cos6» = t aw - r BW (1.1) 

where the r's represent the different surface tensions appearing in the prob- 
lem. This equation can be derived for chemically pure substrates in several 
ways, such as by a mechanical argument relative to the balance of forces, or 
by a thermodynamical argument related to the minimum of the free energy 
of the system ABW @. In these approaches, it is in fact implicitly assumed 
that the surface of the substrate is perfectly flat. If this may well be the case 
at the macroscopic scale (a few mm), it is far from obvious in the presence 
of microroughness. That is to say, how valid this equation will be on top of 
substrates in the presence of atomistic pores or profusions characteristic of 
a solid surface ? 

To examine this question is precisely the aim of our paper. 

Young's equation may be viewed as a direct consequence of the Win- 
terbottom's construction. This construction, first obtained from variational 
principles, describes the equilibrium shape of a crystal as a function of the 
three different tensions that appear in the problem. Its validity at the micro- 
scopic scale, together with the associated contact angle equation, has been 
proved more recently in several models [PD| , |DDR] , PV[| . We consider here 



this construction on top of microscopically rough substrates, in the case of 
1 + 1 -dimensional solid-on-solid models. 

We use two SOS models in fact: one to describe the microscopic interface 
between the substrate W and the fluids A and B, and another one to describe 
the microscopic interface between A and B. For simplicity, we assume here 
that the two models have the same elementary spatial period. 

On the other hand, it is known macroscopically that the roughness of the 
substrate will induce some change in the wall surface tensions, and hence on 
the difference taw — t~bw- This change is described in the literature by the 
so-called Wenzel's equation []W|: 



t aw ~ t bw is proportional to r 
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where r denotes the ratio between the area L of the surface of the substrate 
and that of its projection L on the tangential plane at the contact point 

r = L/L 

We are thus interested to analyze within our model the validity of this 
prediction, extending in that way previous results obtained for the Ising 
model [BDKZ,BDK]. 

In particular, we will also be interested by the corrections to the Wenzel's 
law versus the geometry of the pores or the protusions in our model. This 
research clarifies some preliminary results obtained in that direction with the 
help of numerical simulations [TUBD]. 

The paper is organized as follows. Section 2 is devoted to the presentation 
of the model. Section 3 extends the microscopic validity of the Winterbot- 
tom's construction to rough substrates. In Section 4, we present low temper- 
ature expansions for the wall tension and in Section 5 we compare different 
geometries for our substrate. Concluding remarks are given in section 6. 



2 The model 



To define the model, we consider an SOS model where to each site i of the one 
dimensional lattice we associate an integer variable hi, % — 0, 1, . . . , N, which 
represents the height of the interface between i and i + 1. For a configuration 
h = {ho, ■ ■ ■ , /ijv}j we draw the horizontal lines at height hi between i and 
i + 1 {i = 0, . . . , N — 1), and the vertical lines at each site i, between 



and hi. We use T to denote the corresponding polygonal line (see Fig. [2.1| ). 
Its length is |T| = £ i=1 (l + \ h i ~ h i-i\)- 

We want here to study this interface on top of a rough substrate with 
roughness r. The substrate is thus represented in our case by a periodic SOS 
interface W, with periodicity a, and height configuration h = h , . . . , h^ 
where hi = h n +i, so that 



The energy of a configuration, in a box of length iV (which will be taken 
as a multiple of a), is given by 

H N {r,w) = j AB \r\{rnw)\ + j AW \rnw\ + j BW \w\{rnw)\ (2.1) 
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Here T is above W, which means hi > hi for all i. The set F \ (T n W) 
is relative to the AB microscopic interface, r n W defines the part of the 
substrate in contact with A, and W \ {V fl W) is relative to the contact zone 
between B and W . 

This system describes a system of droplets of a phase B inside a medium 
A on top of the wall W. Jab, Jaw an d Jbw are the energies per unit length 
of the corresponding microscopic interfaces (see Fig. |2.1|) . 




Figure 2.1: A configuration of the interface T on the substrate W. 

Let us first introduce the different tensions appearing in the problem. 

The surface tensions associated to the macroscopic interfaces AB and 
AW are defined as follows : 

t ab (6) = lim log exp(-(3J AB \T\) (2.2) 

N—*oo ply — ' 

where the sum runs over all configurations satisfying h = and = 
iVtan#, and 

r AW = lim log V f exp[-/5F JV (r, W)] (2.3) 

/V — >no /"J /V g ■* 
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where the sum runs over & U configurations such that h = h and h N = 
Tin. Finally, for the interface BW, we have 

T B w = rJ B w (2.4) 

Let us point out that the anisotropy of the SOS model considered here 
leads to an orientation dependent surface tension for the AB interface. That 



the limits exist follows from standard arguments, see e.g. DP, MMR . 



3 Winterbottom's construction on rough 
substrate 

To analyze the microscopic problem of the Winterbottom's construction for 
the model under consideration, we have to consider, following |PDR|| , the 
system given by the Hamiltonian submitted to a canonical constraint 
on the volume of phase B enclosed between the AB and W interfaces. The 
first step consists to analyze the case of a single droplet of the phase B. To 
this end, following |[MK|| , we consider the Gibbs ensemble consisting of the 
configurations which have specified height at extremities and which have a 
specified volume V between the interface AB and the substrate W, 

N 

h N = h N , h = h + M, J2( h i~hi) = V, (3.1) 

8=0 

We assume that the constraints due to the fact that the microscopic interface 
does not touch the substrate 

hi>hi, i = 0,...,N (3.2) 

are satisfied, therefore there is no interaction between the interface and the 
substrate. The corresponding partition function is 

Zi(N, V, M,h)=J2 e^^ (h) S(h N - h N )5(h -h -M) 

U n n (3.3) 

6(J2(hi-h i )-V)Hx(h i >h i ) 

i=0 i=l 

where 



N 

H N (h) = Jab(1 + \hi - h^l) (3.4) 

i=l 
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Hereafter, V and M must be understood as their integer part when they do 
not belong to Z. 

We denote by Z X (N, V, M), the sum of the same Boltzmann factors e~^ HN ^ 
over the configurations satisfying conditions ( |3.1| ) and 

hi>0, i = 0, ...,N (3.5) 



We next consider a conjugate ensemble with the partition function: 

Z 2 (N, u,v) = J2 e ^ HN(h) e W(h)/7V)+ ^ S(h N ) (3.6) 

h 

where V(h) = YliLo ^* an d u E M. and /i6l are the conjugate variables to 
M and V in Z x . 

Our first Theorem establishes the existence of the thermodynamic limit 
for these ensembles and their equivalence in this limit. 



Theorem 3.1 The following limits exist 

^i(u,m) = lim — — \nZx(N,vN 2 ,mN,h) (3.7) 

N^oo [3N 

= lim lnZi(iV,t;iV 2 ,miV) (3.8) 

i/>2(u,n) = lim -— lnZ 2 (N,u,n) (3.9) 

They define the free energies per site associated to the considered ensembles 

as, respectively, convex and concave functions of their variables. Moreover, 
ipi and —ip2 ar e conjugate convex functions: 

—ip2{u, /i) = sup [uv + iim — ipi(v, m)] 

v ' m (3.10) 
i/j x (v , m) = sup [uv + fim + ip2 (u, jj)] 



Proof. We take, for simplicity of the notations, hjsr = and let h mSuX = 
maxj \hi\. Then 

Z X (N, V_, M)e" 2/3|/tmax| < Z X (N, V, M, h) < Z X {N, V+, M)e 2/3| ' imax| (3.11) 
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where 

V ± = V±(N-l)h max 
These inequalities follow by using respectively the changes of variables 

ho — > ho = ho 
hj\i — > hjy = Hn 

hi ^ hi Hi lb /^niax5 ^ 1 5 • • * j N 1 

Now we use the subaddittivity property 

Z X {2N, 2(V + V"), M' + M") > Zt(N, V, M') Z X {N, V", M") e 



(3.12) 

The proof of this property is given in [|MR] and we recall it in the appendix 
for the reader's convenience. From this property, we get the existence of the 
limits when iV — ■> oo of 

(-1//3N) log ZiiN^mN) and (-1//3N) log Z x (N,V + ,mN) 

and the convexity of the corresponding free energy (we take N = 2 n , n G N). 
Provided that h max = o(N) (in fact h max is a constant under the hypothesis 
of section 2) we see from inequality ( |3.11| ) that these limits actually coincide 
with the limit 

ipi(v,m) = lim — — — In Zi(N, v N 2 , mN, h) 
which is thus independent of h, and where 

v = lim — -r = hm — — , m = hm — 

N->oo N 2 N^oo N 2 N^oo N 

We next introduce the partition functions 

Z+(N,u,n)= sup [ e Mvm+W Z X {N,V,M)\ 
v,Mez 

Z*(N,u,n)= [e Pu{v/N)+ ^ M Z^N^M)] 

y,Mez 

and the convex function ^(u, /i) = sup^ m [uv+fj,m—ipi(v, m)\. The Griffiths 
maximum principle adapted to our case (see [M 1\'| |( i.XT^ ) gives that 

lim --^-logZ+(iV,w,/i) = lim -— log Z 2 (N, u, fj) = -^(u,n) 

N-*oo pl\ N^oo pl\ 
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These limits coincide with the limit of (— 1//3N) log Z 2 (N, u, fi), because 
%.{N, u, n)/Z 2 (N, M,/i)^^>0as shown in [PLTRj . QED 

We have thus shown that the free energy of the sessile drop with the 
fixed volume V = vN 2 is proportional to N. Using (|3.10| ), we can determine 



this free energy using the conjugate ensemble where the constraints on the 
volume and on the solid surface do not appear. We shall be interested in the 
case m = which corresponds to a drop on the horizontal plane. We write 

^i0) = ^i{v,m = 0), ip 2 (u) = ip 2 {u,n= -u/2) 

and will express these free energies in terms of the surface tension tab- We 
introduce for that the projected surface tension 

r pr (-tan#) = t ab (6)/ cos 9 

and its Legendre transform 

-(p{x) = sup [xy - r pT (y)] (3.13) 
y 

According to Andreev ||An|| , the Legendre transform ip solves the Wulff vari- 



ational problem when the surface tension is tab{®)- The graph of (p gives the 
boundary of the equilibrium crystal shape of the phase B inside the phase 
A. In the case under consideration, this Legendre transform corresponds to 
the free energy associated with the statistical ensemble conjugate to the en- 
semble defining tab(0) with respect to the constraint on 1%n |MR| (see (|3.23|) 
below). The value of if is given by 

( x 1 , e^ B cosh p J AB 
ip{x) = --log- 



P cosh P Jab — cosh fix 



The following theorem establishes the microscopic validity of the associ- 
ated Wulff 's construction. 



Theorem 3.2 The free energies ip\ and ip 2 can be expressed in terms of the 
functions if and r pr as follows 

u 

i n 

^ 2 {u) = - / tp{x)dx (3.14) 



u 



2 





2 

= — / (p(x)dx — ip(uo/2) (3.15) 



2 



r pT (ip'(x))dx (3.16) 
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where u satisfies 



UQ 

2 1 

<p(x)dx ip(uo/2) = v (3-17) 

Uq J_m Uq 



Proof. Consider the Legendre relation ( 3.1 0|) . The supremum over u, fi is 



obtained for the value w , /io for which the partial derivatives of the right 
hand side are zero: 

v + (d4> 2 /du)(u ,[i ) = 0, m + (d^ 2 /dfi)(u , /x ) = 

That is, for «o,/io which satisfy 

1 f u ° 1 

/ ip(x + fi )dx <p(hq + m ) = v (3.18) 

Uq Jo u 

— [<p(j*q) ~ <p(lM) + wo)] = m (3.19) 

Uq 



The function if is even and thus equation ( j3.19|) for m = gives /x = «o/2; 
inserting this value in ( |3 . 1 8| ) gives expression ( 3.17 ) in the Theorem, and the 
Legendre relation fl3.10| ) reads now 

ipi(v) = uqv + ^ 2 (u ) (3.20) 

if ip'i{v) = Uq or ip' 2 (u ) = —v. The proof now proceeds as the proof of 
Theorem 3 in ||MR| whose main ingredient is the following computation of 
if>i(v). We consider the difference variables 

n i = h i _ 1 -h i , i = l,...,N (3.21) 
and observe that V(h) = ^2f =0 ht = Yli=x ^ n i- Thus 

N 

Z 2 (N, m, AO = II ( S e- /3J - 4s(1+l " l|)+/3(u/7V)ini+/3 ^) (3.22) 

i=l n;gZ 

Since, on the other hand (see |[MR|| ) 

ip(x) = lim — -logVe-^ (h)+ ^5(/i ) = -- In V e -PJA B (i+\n\)+Pxn 
rv ' n^oo 3N ^ v ' 8 ^ 

^ h r nel 

(3.23) 

we have 

N 



Z 2 (N,u,fx) = exp [ - P^2<p(j;i + V) 

i=l 
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and 

N „ N 



Mu, n) = Jim - ^ <p(- t + n) = Jim - ^ -^(-z + /i) (3.24) 

i=l i=l 



y?(x + ji)dx (3.25) 



which for /i = — u/2 gives expression (|3.14|) . The theorem then follows by 
using the Legendre transform relations. QED 



From Theorem |3.2j we know that the equilibrium shape is a piece of 
the associated Wulff shape. To determine which piece, let us introduce the 
ensemble with partition function 

Z 3 (V, At) = e pArN Zi(N, V,M—0) (3.26) 

N 

This ensemble is the conjugate ensemble of the ensemble associated to parti- 
tion function Z\ with respect to the variable N, the conjugate variable being 
the difference of wall tensions At = taw — tbw- 



Theorem 3.3 The limit 

^ 3 (Ar) = Jim -— L= logZ 3 (V, At) (3.27) 

exits and 

^ 3 (Ar) = - sup[aAr - aV>iO~ 2 )] (3.28) 

a 

Moreover 



^s(Ar) = 2 
where Uq is the solution of 



<p(x)dx — uq(p(uq/2) 



(3.29) 



At = ^(u /2) (3.30) 



Proof. The partition function Z3 corresponds to the conjugate Gibbs en- 
semble of the ensemble associated to Z\. Equation ( |3.28| ) is the expression 



of this fact for the corresponding free energies. This can be seen from 

Z 3 (At) ~ Y e ^TN-NM^) „ Bupe-^-Ar^+^^J] 
N 
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which, by taking a = N/y/V, implies equation ( 3.28 ). A rigorous proof 
of equation ( |3.28| ) may be obtained, as in Theorem 3.1, by the Griffiths 
maximum principle (see |GMS||). 



In order to compute ^3 (At) we first remark that from ( |3.28| ) we get 



-a Ar + a ^i(a; c 



(3.31; 



where a is the solution of 

d 

— (-oAt + c#i(a~ 2 )) = -At + ?/>i(or 2 ) - 2a" 2 ^(a" 2 ) = (3.32) 

The function ipx is given by formula ( |3 . 1 5|) which according to equation (|3.17|) 
can be written as 



V>i(f) = 2u v + ip(u /2) 



(3.33) 



where u is the solution of ( |3.17|) . Deriving with respect to v equation (|3.20| ) 
we find 



(3.34) 



Taking a~ 2 = v, equation (|3.32|) reads At + ipiiy) — 2vip' l {y) = which, 
taking ( |3.33| ) and ( |3.34j ) into account, amounts to ( |3.30| ) : At — ip(u /2) = 0. 
Let v be the value of v which according to ( |3.17| ) corresponds to the value 
of u which solves the previous equation, that is 



v = \ [ <p(x)dx - —<p(^) 



(3.35) 



Then from (|3.31|) we get ^(At) = — 4= + -^ip^vo) and, using ( |3.33|) and 
^ 3 (At) = — + -^(2u v + (p(-£)) = 2u \/v^ 



This proves equation ( |3.29| ) stated in the theorem. QED 

Theorem 3.3 gives, by ( |3.30| ), the contact angle relation (see e.g. Re- 
marks 6 and 8 in IIMMKII ) 



r AB (6) cos 6 - t' ab (9) sin 6 = t aw - t bw (3.36) 
This relation reduces itself to Young's equation ( [Lip for isotropic media. 
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Remark 3.1 More general Hamiltonians of the form H = P{\hi — 
where P is a strictly increasing function such that P(x) > \x\, when x — > oo, 
can be treated in the same way as well as the case of continuous height 
variables. 



Remark 3.2 Let us also stress that the proofs may be easily extended to the 
cases of finite range interactions between the interface AB and the substrate 
W. 



Theorem 3.3 establishes the validity of the Winterbottom's construction 
for a single droplet. Actually the initial problem of an interface attracted 
by a substrate with a given volume of the phase B leads to the analysis 



of a gas of droplets with a global volume constraint (see [PDR|| section 4). 



This problem has been considered in ||DDR|| and it has been proven there 
(in the case of a flat substrate) that for a large set of configurations, whose 
probability tends to one when N — > oo, there is only one large droplet, which 
has the macroscopic volume given by the constraint and a gas of microscopic 
droplets without any volume constraint. The free energy of this gas is given 
by the wall surface tension. This shows that the free energy of the large 
droplet is given by Theorem 3.3 and justifies the way we tackle the problem 
provided the same analysis is extended to the case of rough substrate. We 
will not give here the details of such analysis which can be easily adapted 
from section 4 in flBTJR . 



4 Low temperature expansion of the wall ten- 
sion 



This section is devoted to study the behavior, at low temperatures of the 
surface tension taw, defined by equation (|2.3|) . 



Consider a drop of B on the top of the substrate W. Two cases may 
appear: either the liquid B is always in contact with W or there may be 
droplets of A between the liquid and W . Within our SOS model, that means 
that the ground state of the Hamiltonian of the system is given by the mi- 
croscopic interface Y that coincides with the substrate W, or microscopic 
interface T which leave holes between T and W. 
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Let us develop the first case. To this end, we introduce the energy differ- 
ence 

H' N (T, W) = H N (T, W) - H N (W, W) (4.1) 
so that the surface tension taw reads 

taw = rJ A w + hni — -L log Z N (4.2) 

N-*oo pl\ 

where 

Z N = Y,e- pH '» {p > w > (4.3) 
r 

Our first step is to write Z N as the partition function of a gas of elementary 
excitations, simply also called excitations, which can be viewed as micro- 
scopic droplets over the substrate. These excitations are defined as follows. 
Given V and W , we consider the symmetric difference 

a = (r u w) \ (r n w) (4.4) 

We decompose A in maximal connected components A = Si U S 2 U • • • U S n 
called excitations. Here, two components are said connected if they are 
connected considered as subsets of M 2 . A set {Si, 62, ■ ■ ■ , S n } of mutually 
disjoint excitations is called an admissible family of excitations. Then there 
exists a microscopic interface (SOS configuration) T, such that A = Si U S 2 U 
• • • U S n = (r U W) \ (r n W). It is obtained by the formula 

r = (A U W) \ (A n W) (4.5) 

This correspondence between admissible family of excitations and interfaces 
SOS configurations is one-to-one. 

The energy difference H' N reads in terms of families of excitations as 
H' N (T,W)=E(Si) + --- + E(S n ) 

where 

E(S) = J AB \S\ (S n W)\ - {Jaw ~ JbwW n W)\ (4.6) 

Indeed 

K = j ab \t \ (r n w)\ + j aw \t nw\ + j bw \w \ (r n w)\ - J AW \w\ 

= J AB \T \ (r n W)\ - {Jaw - JbwW \{mw)\ 
= J AB \A \ (A n W)\ - {Jaw ~ Jbw)\& n W\ 
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The equality r\(rnW) = A\(AlW) follows from equations Q, (gj^) and 
TUW = AUW. The equality |W\(rnW)| = |AnW| follows from equation 
JOP and the relations r U W = A U W, |W| + |T| = j W U T| + \W D T|. 

Then 

n 

Z* = £ JJ e -/H5(*) (4.7) 

A={<5i,...,5 n }cA ]v i=l 

where the sum runs over admissible families of excitations whose projection 
is included in the infinite cylinder A^ = {x = {x\,X2 ) e 1R 2 : < xi < N} 
and the product is taken equal to 1 if A = 0. 

In the concept of excitation that we are considering, the configuration 
T = W, in which the microscopic interface is following the wall, as the 
ground state of the system. In other words, we assume that H' N (T, W) > 
for all T and N, or equivalently, that 

mmE(5)>0 (4.8) 

5 

In fact it is enough that this condition is satisfied for N = a, that is for all 
excitations belonging to A a . 

We next consider arbitrary families of elementary excitations non neces- 
sarily mutually compatible and in which a given excitation can appear several 
times. To any such family {Si, . . . , S n } a graph G(di, . . . , S n ) is associated in 
such a way that to each excitation corresponds (in a one-to-one way) a vertex 
of the graph, and there is an edge joining the vertex corresponding to Si and 
Sj whenever Si and Sj are not compatible or coincide. We introduce the clus- 
ters C as the arbitrary families of excitations for which the associated graph 
G(di, . . . , S n ) is connected (this means that the excitations draw a connected 
set in H 2 ). Then we get 

\ogZ N = J2$ T (C) (4.9) 

c 

where the sum runs over all clusters whose excitations belong to A^. The 
truncated functions $ T are defined by 

® T (Si, . . • , S n ) = a( ^'V A) n z~ mk) (4-10) 



the arithmetic coefficient being 
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a(8i, . . .,8 n ) 



GcG{S U - 



:-i) 



1(G) 



(4.H) 



,5n) 



Here the sum runs over all connected subgraphs G of G(S\, . . . ,S n ), whose 
vertex coincide with the vertex of G(5\, . . . , 5 n ), and £(G) is the number of 
edges of the graph G. If the cluster C contains only one excitation then 
a(S) = 1. 



To express condition Q4.8Q in terms of the coupling constants, we need a 
description of the substrate. Let T(z) be the horizontal line at height z, that 
is hi = z for all i. For any z G Z such that infj ^ + 1 < z < sup^ the 
substrate IV and the line T(z — e), < e < 1, intersect in a finite number of 
points, PI r(2 — e) = {At, A 2 , . . . , A p }, ordered in such a way that the first 
coordinates i k (k — 1, . . . ,p) of Afc satisfy z'i < %i < ■ ■ ■ < i p . The part of W 
between the two points B k = (ik,z) and B k+ i = (ik+i,z) lies either below 
either above or on the the substrate W. It is called a well in the first case and 
we denote it by w k (z), and a protusion in the second case (see Fig. We 
let p = max^fc \w k (z)\/(i k+1 - i k ) = max Zjfc \S k (z) H W\/\6 k (z) \ (6 k (z) H W)|, 
where is the excitation 8 k (z) = w k (z) U [i k ,i k+ i\. 



W 



_ l_Bi 

r(z-e) 



5, 



5 4 



» ! 3 



«2 



*3 



Figure 4.2: The wells wi between B\ and £?2 and W3 between B 3 and £> 4 . 
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Then condition 



reads 



Jab > p{Jaw — Jbw) 



(4.12) 



Hereafter it will be more convenient to denote W the infinite periodic wall 
whose restriction to An is given by the previous height ho, ■ 



. ,/ijv- Notice 

that the expression ( [4.6D of the energy of excitation remains unchanged. We 
shall use W a to denote the restriction of W to A a . 



Theorem 4.1 Assume that the condition ( 4.12| ) is satisfied, then, for any 
(3 > p = 1.9 (1 + p)[Jab — p(Jaw ~ Jbw)] 1 , the following serie, defining 
the wall-medium surface tension, is absolutely convergent 



TAW = rJ AW --^^ 



beWa CBb 



$ r (C) 

icrwi 



(4.13) 



Proof. The proof of formula ( J4.13 ) as well as that of the absolute conver- 
gence of the series can be established following ||GMM|| (Chapter 4) in which 
the low temperature contours of the Ising model were considered in the role 
played here by the excitations (see also |D|, [KP| , 



The first ingredient is the following lower bound on the energy: 

E(5) > {1 + P )- 1 [Jab-p{Jaw-JbwW\ 
This bound follows from expression 



(4.14) 



E{8) 



J 



Jab + Jaw — Jj 



BW 



AB 



1 + 



|<5\(<5fW)| 
\5r\w\ 



obtained from ( glj ) and the inequality (|5rW|/|<5\(<5rW)|) < p that holds 
true for any excitation 5: this inequality is a consequence of easy geometrical 
arguments used with condition ( |4.12| ). 

Inequality (|4.14j) together with the fact that the number of polygons (or of 
excitations 5) of length I passing to a given point is less then 3 e ensures in 
particular the convergence of the series Y2sBb e~ l3E ^ , for all bond b as soon 
as (3 equals some (3' Q . 

The convergence of the cluster expansion needs furthermore the existence of 
a positive real- valued function fj,(5) such that 



-^Vr X exp{ ]>>(<*)} 

^ S'iS ' 



< e~ a < 1 



(4.15) 
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where the sum runs over excitations 5' incompatible with 5: this relation is 
denoted by 5' i 5 and means that 5' do not intersect 5. Taking in addition with 
the above remark on the entropy of excitations, that the lengths of excitations 
are even with minimal value |5 min | = 4, that /-t(<5) < \S\ X^'gbM^'X an d 
choosing p(5) = (Se*) - ' 5 ', inequality ( |4.15| ) will be satisfied whenever 



p (1 + p)- x [ J AB - p(J AW - J BW )} > log 3 + t + 



-it 



1 - e 



~2t 



The value to ~ .61 that minimizes the function t + [e 4t / (1 — e 2t )} provides 
the corresponding (3q given in the theorem. The expression (gjj) then follows 



from (|4.2|) and Q4.9|) by setting N = pa and letting p — > oo. QED 

The other cases where the ground state is not the wall W will be discussed 
elsewhere |[DMR| . They lead in particular to the Cassie's law Q. 



5 Comparison of different geometries 



We shall now restrict to some specific walls W. Namely, we assume hi = 
for % = c — 1 and hi = b for i = c a — 1, (1 < c< a — 1), see 



Fig. p 



We will denote by Ar(r, c) = — r BH/ the difference between the 
surface tensions corresponding to the roughness r and the parameter c. The 
roughness has value r = 1 + 2b /a and it is independent of c. 

Let us stress that Theorem [4.1| implies Wenzel's law at low enough tem- 
perature 

Ar(r, c) = t-(Jaw — Jbw) + corrections 

The next theorem concerns the corrections to Wenzel law. It gives a 
comparison between different geometries with the same roughness by varying 
the parameter c. 

Theorem 5.1 Assume that J = Jab > and J' = Jaw ~ Jbw > satisfy 
J - (2b + I) J' = 2(6 + l)K x > 0, then, for (3 large enough, f3M > {b + 
l)(1.9(a + c) + .56) where M = min{(6 + 1) J', 2(b + 1)7^, | J' - (b + 



a) if 2 < c < a - I, 



Ar(r, 1) < Ar(r, c) 
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b) */l<c<Co-l, 



c) if Cq < c < a — 2, 



Ar(r, c) < Ar(r, c + 1) 



Ar(r, c) > Ar(r, c + 1) 



where when a is odd Cq = (a + 5)/2 and when a is even cq = a/2 + 2 if 
J' <{b + l)Ki and c = a/2 + 3 if J' > (b + 1)^ . 



This result is represented graphically in Fig. 5.1 




Figure 5.1: Plot of the wall tensions difference Ar(r, c) as function of the parameter c. 



It means that for a given roughness r, there is an optimum in the wall tension 
at c = Co- That is to say that for c = Co, the associated contact angle 9 for the 
sessile drop will be minimum for 9 < tt/2 and maximum for 9 > n/2. These 
results also confirm the data obtained by numerical simulations in ||TUBD 



It was indeed observed that Ar(r) for 2c? random substrates with a fixed 
roughness r remains between the "single" protusion and "single" hole case 
as represented in Fig. f>72. 
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0.7 - 



0.6 1 1 1 1 1 1 1 1 1 

0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4 



Figure 5.2: Plot of Ar(r) as function of the roughness r at two different temperatures, 
for the single hole c = 1 (lower dashed lines), the single profusion c — a — 1 (upper 



dashed lines) and the random geometry (diamond), [TUBE] 



On the basis of Theorem |5.1| , we have in fact 

Ar(r, c = 1) < Ar(r, c) < Ar(r, c ) 
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Since on the other hand we have 

Ar(r, c = a — 1) ~ Ar(r, c = cq) 

we can understand that indeed single protusions and single holes will be good 
approximations for the upper and lower limits of the wall tension as already 



indicated in TUBD 



Let us also point out here that the numerical simulations seem to indicate 
that the first order correction to Wenzel is already enough to describe the 
wall tensions up to the half of the 2d Ising critical temperature. 

That this optimal geometry also holds for more general systems remains 
up to now an interesting open question. 

Before proving Theorem we give the following 



Lemma 5.1 Assume that cJ — (c + 2b) J' = 2(c + b)K c > 0, then, for f3 > 
2(1.9 + a/4) (b + c)[cJ- ( y c + 2bJ')]-\ a > we have, 



$ T (C) = ll + e(c)]e 



■P[cJ-(c+2b)J'] 



(5.1) 



C3(0,0) 
C3(0,c) 



where 
\e(c)\ < e 2u{c+b)+ ^ 



1 _ e -/3(J+J') 



-P(J+J') 



-2(J3J'+u) 



+ 



1 _ e -2(pj'+u) 



v = 1.9 + a/4, 1] = log 



-4t 



1 - e-*o 1 - e- a 



log 



0.2 e~ 

1-e- 



Proof. We first observe that all excitations satisfy maxs(\5 C\W\/\5\ (5 C\ 
W)\) < (c + 2b) /c, i.e., p = 1 + 2b/ c, so that inequality flQ3|) reads 

Under the condition, /3X C > 1.9+a/4, a > 0, the cluster expansion converges 
and moreover 

\^ T (C)\ < exp{(-/3K c + u)m + 7]} (5.3) 

\C\>m 
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where for a cluster C = {5i, . . . , 6 n } we use |C| = \Sx\ + • • • + \S n \ to denote 
its length. To prove (|5.3|), we write for a cluster C = {81, . . . , 5 n } of length 
at least m, |$ T (C)| < {l/n\)e-^-^ m a{5 1 , . . . , 5 n ) Uti e~^ E ^\ Then 
we use that condition ( |4. 15| ) is satisfied if for any 5, /3K C > log 3 + t + 
^_ e -2t + ot/\8\, i-e. for [3K C > 1.9 + a/4 by choosing /i and t as in the proof 
of Theorem [O]. Under ( ^.14| ), one knows, c.f. pMM|| , that £ C9(5 $ T (C) < 
/i(5)e- Q /(l - e-°). We use finally £ feb //(<*) < e" 4 7(l - e" 2 *) = .2 when 
t = .61 to get flO|). 

We let 5o be the excitation corresponding to the interface r , 5q = (r U 
W) \ (r U W), where T is given by the height hi — b for i — 0, 1, . . . , c — 1 
and /ij = ft,, otherwise. That is <5o is the boundary of the rectangle R = 
{(x,y) 6l 2 : < x < c, 0, < y < b}, see Fig. |1| Its energy is: E(5 ) = 
cJ — (c + 26) J'. Denote by W(bi, b 2 ) the part of the wall between the points 
B 1 = (OA) and B 2 = (c,b 2 ). Then, 

<$>T {c)=e -picj- {c+ 2 b )j> ]+ $ T (C) (5.4) 

C:CnWDW(b,b) C:CnWDW(b,b) 

\C\>2(c+b+l) 

Indeed the first term of the R.H.S. of Q5.4p corresponds to the excitation 
5o and the second terms run over the other clusters containing W(b, b) the 
length of them being at least \8q\ + 2 = 2(c + b + 1). By (|5.3|), this term is 
bounded as follows: 

\^ T (C)\< exp{-2(c + b+l)(pK c -u)+r]} (5.5) 

C:CnWDW(b,b) 
\C\>2(c+b+l) 



Next we observe that for all excitations whose intersection with the wall 



is W(b u b 2 ) satisfy m.ax i {\5nW\/\8\{5nW)\) < (c + b x + 6 2 )/(c + |6i - b 
Hence the bound (|5.2|) is improved as follows: 



2 ■ 



E(6) > (C+ ' 6l - fe2l l J 'i,;\ + ,? 2)J V l (5.6) 
(2c + bi + b 2 + I&! - 6 2 |) 

when i5nF = 6 2 )- All the associated clusters have length \C\ > 

2c + b\ + b 2 + I b\ — 62 1- Therefore 

^ |$ T (C)| < exp{-/3[(c + \b x - b 2 \) J - (c + h + b 2 )J'}} 

C:CnW=W(b u b 2 ) 

x exp{(2c + 61 + b 2 + \b x - b 2 \)v + 77} (5.7) 
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Thus, 



J2 Yl |$ T (C)|<2[l-e-W)] 



M-l-l 



6l+6 2 <2b-l C:CnW=W(b 1 ,b 2 ) 
6i>l,b 2 >l 

x (exp{-/?(c + l)J + P(c + 2b- l)J' + 2v(c + b) + n} 

+ [1 - e - 2 ^ J>+ ^] ~ X exp {-(3cJ + (c + 2b - 2) J' + 2z/(c + 6 - 1) + r?} 

(5.8) 

The first term inside the parenthesis comes from the summation over 1 < 
6 2 < b\ — 1, bi — b and the second term from the summation over 1 < 6 2 < bi, 
1 < h < b- 1. 

Using that 

C9(0,0) C:CTWDW(b,b) 6i+fe 2 <26-l C:CnTy=U / (6i,6 2 ) 

C9(0,c) bi>l,fe 2 >l 



the proof follows from ( |5.4|5.5jj5T8"|) . QED 

Proof of Theorem [37 



The lower bounds ( |5.2y5.6|) on the energy can be improved for some excita- 
tions. Let A(i, j) = {x = (x, y) G 1R 2 : i < x < j} denote the infinite cylinder 
between the vertical lines x = i and x = j. For the excitations included in the 
strips A(l -a, c-1) and A(l, a+c-1), one has max s (\5nW\/\5\(5nW)\) < 1, 
and thus by arguing as in the proof of ( |4.14| ) 

E{5)>\{J-J')\5\ (5.9) 
The associated clusters satisfy thus: 

\® T ( C )\ ^ ex p{ ~ f( J " J')m + vm + v} (5-10) 

\C\>m 
C3b 

For the excitations included in the strips A(l,c — 1) and A(c, a), one has 
\5nW\ < \S\/2- 1. Therefore 

E(S) = J\S\ - (J + J')\S HW\> -(J - J')\S\ + J + J' (5.11) 
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and the associated clusters satisfy: 

\® T ( C )\ < ex P { - ~( J -J')m- f3(J + J') + urn + r/ j (5.12) 

\C\>m 
CBb 

We let 5\ be the excitation corresponding to the interface Tx, S% — (I\ U W) \ 
(Ti U W), where T\ is given by the height hi — b + 1 for i — c — 1, . . . , a — 1, 
and hi = hi otherwise (Fig. |5.3|) . Its energy is E(Si) = (a — c + b + 3) J — 
(a-c + b + I) J' and \5x\ = 2(a + b + c + 2). 

We let S^ k \x,y) denote the excitations of width k, height 1 (and length 
2(k + 1)) whose intersection with the wall is the segment [(x,y), (x + k,y)] 
(of length k). Their energy when y = b or when y = and they do not 
intersect the vertical part of the wall are: = k(J — J') + 2J. 




Figure 5.3: The excitations <^ fc \ and Si translated by —a. 



Then we the decompose the sum involved in (|4.13j) as follows: 



£ $ T (C) =S 1 (c) + 5 2 (c) + 5 3 (c) 

CnA(O,a)^0 



(5.13) 
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where 

CcA(l,c-l) CcA(c,a) 
CcA(l-a,c-l) CcA(l,a+c-l) 

cn(o,o)^0 cn(o,c)^0 

5s(c)= £ ^ 

cn(o,o)^0 
cn(c,o)^0 

Let us compare the differences Si(c) — 5j(c + 1), i = 1, 2, 3. When a is even, 
we have 



S 1 (c)-S 1 (c+1) 



where 



e -P<J~J')-WJ + R X ( C ) if c _ 2 < a - c - 2 

_ e -/9(a- C +l)( J- J')-2/3J + i^( C ) if C - 2 > a - C 

(5.14) 



|i?l(c)| < 4 £ < 4e -/3( C +2)(J-J')-/3(J + J')+2(c + 2),+r / 

CcA(c,a) 
|C|>2(c+2) 

|i?i(c)|<4 V] < 4e -/3(a-c+3)(J-J')-/3(J+J')+2(a-c+3)^+r, 

CcA(l,c) 
|C|>2(a-c+3) 

(5.16) 

Indeed, there is a one-to-one correspondence between the clusters C of base 
of size |C fl W| = occurring in «Si(c) and Si(c + 1) till reach some value. 
This value, when c — 2 < a — c — 2 is precisely c, because in that case there are 
clusters (of base of size c) which belong to A(c, a) hot neither to A(l, c — 1) 
nor to A(c + 1, a). There is precisely one excitation 8^ c ' of base of size c (and 
length 2(c + 1)) which belong to A(c, a) bot neither to A(l,c — 1) nor to 
A(c + 1, a). Its energy is = c(J — J') + 2J and gives the corresponding 
term in ( |5.14j ). The other clusters have length |C| > 2(c + 2). This gives the 



first bound on the reminder i?i(c). The second bound in ( |5.15| ) follows from 
( |5.12| ). When c — 2 > a — c the argument works in the opposite direction. 



The value k is a — c + 1 and there is a corresponding <jw (of length 2(a — c+2) 
which belong to A(l, c) but nor to A(l, c — 1) nor to A(c, c — a). Its energy is 
(a — c + 1)( J — J') + 2 J and provides the corresponding term in ( |5.14j ). The 



other clusters have length \C\ > 2 (a — c + 3). This gives the bound on the 
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reminder R[{c), the second inequality in Q5.16| ) following from ( p,12| ). 
When a is odd: 

( e -f3c(J-J>)-2(u + Ri ( c ) if c _ 2 < a - c - 3 

Si{c) - Sx(c+ 1) = < ifc-2 = o-c-l 

_ e -/3( a - c +l)( J- J')-2/9J + ( c ) if c _ 2 > a - C + 1 

(5.17) 

where must be understood as YliOM & T (C) with M as large as we wish. 
Indeed the same reasoning as for a even applies. In addition there is the 
particular case c — 2 = a — c — 1 where the width of the cylinder A(l, c — 1) 
equals the one of A(c — a — 1) and the width of A(c — a) equals the width of 
A(l,c). 

For ^(c), we have: 

\S 2 (c)-S 2 (c+l)\ = \R 2 (c)\<2 l$ T (C)l+2 Yl 1^)1 

cn(o,c)^0 cn(o,c)^0 

|C|>2(c+l) |C|>2(a-c+6+2) 

< 2e~ l3( - c+1 ^ J ~ J '^ l+2( - c+1 ^ u+v + 2e~ /3l - a ~ c+b+2 ^ J ~ jn, ~ l3( - J+J ^ +2( - a ~ c+b+ ' 2 ^ u+ ' n 

(5.18) 

Indeed the clusters of minimal energy containing (0, c) and for which the 
correspondence is not one-to-one are the excitations S^ c \0, 0) of length 2(c + 
1) and the excitation 5i of length 2 (a — c + b + 2). All other clusters have 
length greater or equal than either |5^(0, 0)| + 2 or than |d>i| + 2. To bound 
the first sum we used ( p-l(J| ) an d to bound the second sum we used ( |5.12| ). 



Finally, by Lemma |5.1| 

S 3 (c) - S 3 (c + 1) = (1 + £ >) e -PlcJ-(c+2b)J>] ( 5 19) 

where 

\e'\ < |e(c)| + e-^ J ~ J/) (l + |e(c + 1)|) (5.20) 

Assume now that 2c > a + 2 if a is even or 2c > a + 3 if a is odd, then 
from (|5TTi )-(|CTD 



Ar(c + 1) - Ar(c) = ^ Si(c) -Si{c+ 1) 



i=l 



1 + £i)e -/3(a- C+ l)(J-J')-2/3J + (1 + £2)e -0c(J-J>)+2pbJ> ( g 21) 
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where 



| £l | < 4 e -^(J-^')+2(a-c+3)i/+7 ) + 2 e -Pb(J~J')+2(a-c+b+2)u+ri (5.22) 

\e 2 \ < \e'\ + 2e-^- J ')-2/3W+2( c+ 2),-H ? (5 _ 23 ) 



Therefore for (3 large the sign of ( 5.21 ) will be given by the sign of the 
difference 

E (a-c+i) _ E ^ = ( a _ c + X )(J _ J/) + 2 j _ C (J _ J') + 2bJ' 
More precisely 



Ar(c+ 1) - Ar(c) < A 
a + 5 



a + 4 



— c — 



1+M 



.4 



2 ./ - ,/' 2/3(J - J') 



— c — 



1)^ + 2(6 + 1)1^ !ogT^g 



J - J' 



2(3{J-J>) 



At(c+ 1) - Ar(c) > 5 



a + 4 J'- (6+1)*! ^SiqS 



j- j' 



2(3(J-J')_ 



q + 3 (6 + l)J' lo §™ 



J - J' 2/3( J - J') 



(5.24) 
(5.25) 
(5.26) 
(5.27) 



where A = 2(3(J - J')(l + | £2 |) e -^( J - J ')+2/3W and fi = 2 /?(j - J')(l + 
|£i|)e~ /3 ( a ~ c+1 * )< '' 7 ~' 7 )~ 2/3J . On the other hand when 2c < a if a is even or 
2c < a + 1 if a is odd, we get from (|5.14|H|O0D 



Ar(c + 1) - Ar(c) > e~^ J - J>) 



JPbJ', 



l-|e , |) + e- 2/3J x(2c<a+l) 



(5.28) 



Therefore, for a odd, inequality (|5.25|) proves Statement c) of the theorem 
while (|5T27|) gives Ar(^) < Ar(^) i.e. Statement b) for c = (o + 3)/2 the 
result for the other value following from ( 5.28|) , all this provided /3 is large 
enough as stated in the hypotheses of the theorem. When a is even we can 
conclude only if J' ^ (b + In this case, the statement c) follows from 

( 5.24 ) while the statement b) follows from (|5.26|) and (|5.28 ). 



To find the lower bound on (3, we let e\((3M) and £%(f3M) be respectively 
the upper bounds on \ei\ and |e 2 | obtained by replacing J' by M/(fe + 1) 
and K x by M/2{b + 1) in fl5^2]) and (£2^). We use also K c - vc > M/(b + 
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1) - (a - \)v. Then, we take f3M/(b + 1) = (a + b)v(o) + (l/2)(r/(a) + a). 
The value of a giving the lower bound stated in the theorem ensures that 
(l/2/3M){log[l+£i(/3M)]-log[l-£ 2 (/5M)]} < 1, e x ((3M) < landei^M) < 
1. 

The above analysis leads also to: 



-f3(J-(l+2b)J') 

Ar(r, 1) = r(J AW - J BW ) — + O^"^ 26 -^ >) (5.29) 

Ar(r, a - 1) = r( J AW - J w ) - X {a > 3) + O^ 3 ^" J >) 

pa 



(5.30) 



The second term in the R.H.S. of ( |5.29| ) comes from the energy of the ex- 
citation 5 for c = 1 and the second term in the R.H.S. of ( |5.30| ) comes 
from the energy of the excitation 5^^(0,0). Relations ( |5.29| ) and ( |5.30| ) give 
Statement a) and end the proof of Theorem |5.1| . QED 



Notice that by (|5.29| ) the first order term for the corrections of Wenzel's 
law is given in the case c = 1 by (— 1/ [3a) exp{/3[a(r — 1)(Jaw — Jbw) — Jab + 
Jaw — Jbw]} and thus decreases with the roughness r. 



6 Concluding remarks 



Within the 1 + 1-dimensional SOS model, we have shown that the shape of 
a sessile drop on a rough substrate is given by the Winterbottom's construc- 
tion. Using low temperature expansions we have analyzed the first order 
corrections to the Wenzel's law describing the contact angle 6 of this ses- 
sile drop. These corrections vanish as the temperature goes to zero. For a 
given family of two steps substrate (0, 6) with the same roughness, we have 
shown that there is an optimal geometry that maximize or minimize the wall 
tension. These results confirm and explain previous numerical simulations. 
That is to say that, at least for this given family of substrates, there is indeed 
an optimal geometry for wetting. That this property can be proved for more 
realistic models or experimentally is an open challenging question. 
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Appendix 

Lemma The partition function Z\ satisfies the subadditivity property 

Zi(2N, 2(V + V"), M' + M") > Z X {N, V , M') Z 3 (N, V", M") e" 2/3 ^ 

(A.l) 



Proof. In order to prove this property we associate a configuration h of the 
first system in the box of length 2N ', to a pair of configurations h' and h" of 
the system in a box of length N, as follows 

h 2i = ti i + h", i = 0,...,N . 
h^^hl. + h't ,t = l,...,N 

Then h 2N = h' N + h" N = 0, h = h' Q + h% = M' + M" and 

N N N 

i=0 i=o i=l 

= 2 [V(ti) + V(h")} - M" 

This shows that the configuration h belongs to Z X (2N, 2{V'+V")-M", M' + 
M"). Since H 2 ^{b) = -Hjv(h') + Hjv(h"), because n 2i = n\ and n 2 i-i = n" 
for i — 1 , .. . , N — 1, as follows from ( |A.2| ) , we get 

Z X {N, V, M') Z 3 (N, V", M") < Z X {2N, 2(V + V") - M", M' + M") 

Then we use the change of variables hi = hi + [M"/(2N — 1)] for % = 
1, . . . , 2N — 1, ho — ho, h 2 N = h 2 N — which gives 

Z X (2N, V - M", M) < e mM"\/{2N-i) Zi{2N ^ M) 
to conclude the proof. QED 
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